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We report on a theoretical study of the electronic, optical and magneto-optical properties of dig¬ 
ital ferromagnetic hetero-structures based on Mn (5-doped GaAs. We consider different structures 
corresponding to Mn contents within the range 12-50% and we study how the system changes as a 
function of the doping concentration. Our first-principles approach includes the spin-orbit interac¬ 
tion in a fully relativistic pseudopotential scheme and the local fields effect in the description of the 
optical absorption. We show that Mn (5-doped GaAs shares many properties with the uniformly 
doped Gai-xMnxAs system, i.e. half-metallicity, similar absorption spectra, and moderate Kerr 
rotation angles in the visible spectral region. 
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I. INTRODUCTION 

Digital Ferromagnetic Heterostructures (DFH) have 
been recently proposed as a way to overcome some 
of the drawbacks of random-doped Gai-xMuxAs Di¬ 
luted Magnetic Semiconductors (DMS), without loos¬ 
ing the peculiar properties of the latter connected with 
the possibility to exploit spin-polarization effects within 
a semiconducting (or half-metallic) medium compatible 
with the present technology of semiconductor electronic 
devices^. Mn-based DFH in HI-V semiconductors are 
super-lattice heterostructures where (fully or partially 
filled) mono-layers (ML) of Mn -which is substitutional 
to cation group HI elements- are spaced by a fixed num¬ 
ber of layers of the HI-V compound. On one hand, these 
(5-doped structures allow a concentration of Mn higher 
than in random-doped (Mn-)HI-V semiconductors, and 
thus, due to the different valence between the group HI 
atom and the substitutional Mn, can produce in principle 
a larger concentration of holes, enhancing the Curie tem¬ 
perature (Tc) of the system. On the other hand, DMS Tc 
can increase considerably also because it depends on the 
spatial distribution of magnetic ions, which can change 
substantially in low dimensional systems: Nazmul and 
CO- workers^ showed that in a p-type selectively doped HI- 
V hetero-structured composed of Mn (5-doped GaAsjp- 
AlGaAs a remarkably high Tc of about 250 K can be 
reached; in Ref. Q Chen et al. report Tc well above 
room temperature in GaSb/Mn digital alloys. 

Possible applications of DMS range from spintronics 
(i.e., the idea to use the spin degree of freedom in or¬ 
der to transmit, write, or store digital data)^i^ to the use 
of spin-polarized currents for the study of surface mag¬ 
netism and magnetic anisotropies in materials. In this 
respect, half metals, which are conductive in one spin 
channel and semiconducting in the other, are ideal can¬ 
didates to build spin-injectors. Many DMS have indeed 
an half-metallic behavior—. Moreover, several DMS also 


display intriguing magneto-optical coupling effects: an 
example is the Magneto Optical Kerr Effect^ (MOKE), 
i. e. the polarization rotation of the light reflected by 
the surface of a magnetic material, which is at the basis 
of the operation of magneto-optical data storage media 

Moreover, the MOKE can be used in the study of elec¬ 
tronic structure in magnetic materials. 

Previous studies of GaAs based DEH include 
the experimental realization, by means of low- 
temperature Molecular Beam Epitaxy (MBE), of 
GaAs I Gai-xAsMnx, separated by GaAs spacers of vari¬ 
able thicknesses^, and measurements of their optical con¬ 
ductivity spectra in the infrared regioit^. 

On the theoretical point of view the optical and 
magneto-optical properties of Mn-doped GaAs have 
been studied by means of the linearized Eull-Potential 
methods^iiS (ELAPW approach) or by the Korringa- 
Kohn-Rostoker methodii, but only in the case of an 
uniform Gai-xMnxAs doping. Ab-initio theoretical pre¬ 
dictions for the magneto-optical properties of DEH are 
hence lacking. 

A possible advantage of DEH J-doped structures with 
respect to uniformly doped Gai-xMnxAs DMS derives 
from the possible amplification of MOKE effects by ex- 
citon confinement. However, also local field effects (cor¬ 
responding indeed to the electron-hole exchange in the 
excitonic Hamiltonianpi^ ) may play an important role 
in such a strongly anisotropic system. However, the 
intrinsic complexity of DEH, with unit cells contain¬ 
ing many atoms, makes the direct application of full- 
potential methods such as those of Refs. and com¬ 
putationally prohibitive. 

On the other hand, more efficient ab-initio methods 
based on pseudo-potentials (PP) and plane waves have 
not been systematically used for the study of MOKE 
spectra, because of the common wisdom that an all¬ 
electron approach was required in order to describe the 
wave-function within the core region (where the spin- 
orbit (SO) interaction is stronger). Instead, some of us 
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Figure 1. (color online) Atomic structure (unit cells) of some 
of the studied (5-doped materials. Yellow spheres are As 
atoms, the violet ones represent Ga, and the cyan color is 
used for Mn. The corresponding Mn-doping concentrations 
are: 12.5% (rightmost structure, an 8-layers super-cell with 
7 GaAs layers separating MnAs ones); 16.7% (6-layers super¬ 
cell with 5 GaAs layers); 25% (leftmost, 4-layers super-cell 
with 3 GaAs layers); 50% (2-layers super-cell, not shown). 


Figure 2. (color online) Spin-resolved Density of states of 
the four studied Mn-doped GaAs heterostructures, compared 
with that of bulk GaAs (black dashed line). Panel a to d cor¬ 
respond to decreasing Mn concentrations, from 50% to 12.5%. 
Spin up and spin down components are computed by project¬ 
ing the spinorial wavefunctions onto the S^ spin eigenstates. 
The inset of each figure shows a zoom in the vertical scale of 
the density of states in the minority spin channel close to the 
Fermi level. 


have recently demonstrated that, provided that the spin- 
orbit interaction is fully taken into account in the con¬ 
struction and use of “fully relativistic” PP, the plane- 
waves based approach provides results of the same qual¬ 
ity as those of all-electron calculationsi^. 

Indeed the MOKE is due to a different dielectric re¬ 
sponse for right (e+) and left (e_) circularly polarized 
light, which stems from the difference between the wave- 
functions with -\-Lz character and —Lz character. In a 
magnetic material this difference is due to the SO cou¬ 
pling in the Hamiltonian, H = • S, which transfer the 

spontaneous symmetry breaking between the -\-Sz and 
the —Sz spin component to the spatial part of the wave- 
function. Thus the MOKE cannot be captured by a per¬ 
turbative description of the SO effect, nor by the use 
of scalar relativistic PP, which only include the effect of 
the SO in the energy levels without affecting the wave- 
function. 

In this work we hence perform first-principles simu¬ 
lation of Mn (5-doped GaAs structures^, including the 
effect of the SO in a non perturbative way and using 
fully-relativistic PP. We then compute the absorption 
coefficients of the system and the MOKE parameters in 
the so-called polar geometry, where the magnetization 
vector is oriented perpendicular to the reflective surface 
and parallel to the plane of incidence. This is the most 
studied geometry and the one which gives the largest 
Kerr rotation. 

This work is organized as follows. In Sec. [Ill we present 
DET-LDA ground state results and Kohn-Sham (KS) 
band structure for heterostructures described by 8,6,4, 
and 2-layers super-cells (SC), corresponding to doping 


concentrations from 12.5% to 50% (see Eig.[T]). We show 
the electronic density of states, projected onto spin com¬ 
ponents, atomic species, and s/p/d angular momenta. In 
section im we deal with optical absorption spectra, com¬ 
puted within the linear response scheme in the random- 
phase approximation (RPA), and compared with experi¬ 
ments. We then extend our analysis to magneto-optical 
effects in section lTlI Ai after briefly summarizing the basic 
equations for the MOKE and MCD spectra, and we dis¬ 
cuss our results analyzing the similarities and differences 
with the case of uniform doping. Einally, in section IIVI 
we draw our conclusions. 


II. ELECTRONIC PROPERTIES OF Mn DOPED 
GaAs: DELTA DOPING 

We model Mn (5-doped GaAs by means of the periodic 
super-structures shown in Fig. [TJ The variable thickness 
of bare GaAs interposed between MnAs layers allows us 
to simulate (5-doped materials with 12%, 16%, 25% and 
50% of Mn atoms replacing Ga ones, (with the latter case 
corresponding to an hypothetical GaAsjMnAs super¬ 
lattice) . 

We perform our calculations with DFT-LDA-i^“— in 
the plane waves PP scheme, as implemented in the Quan¬ 
tum Espresso (QE) packagei^. To keep into account the 
spin-orbit effects on the valence wave-functions, the use 
of fully relativistic (norm conserving) PP is necessary. 
Moreover, the presence of non-collinear spin polariza¬ 
tion (i.e., spinorial KS wavefunctions) must be allowed. 
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Figure 3. (color online) Kohn-Sham band structure and total 
DOS for the 25% doping case (panel (b) of figure [2|). The red 
line is drawn at the Fermi energy. 


Figure 4. (color online) Kohn-Sham band structure and total 
DOS for the 50% doping case (panel (a) of figure [2|). The red 
line is drawn at the Fermi energy. 


Both features are available in QE. We use tetragonal unit 
cells, by relaxing the lattice parameter along the growth 
direction (perpendicular to the Mn doped layers) and all 
the atomic positions within the unit celU^. The in-plane 
lattice parameter is kept fixed at its bulk GaAs value, 
consistently with an epitaxial growth on the GaAs (001) 
oriented substrate. Ground state calculations are per¬ 
formed with a kinetic energy cut—off of 50 Ry (680 eV) 
for the wave-functions, and a Monkhorst-Pack^^ grid 
12xl2xn for the Brillouin zone (BZ) with n = 8,4, 3, 2 for 
the two, four, six and eight layers SC respectively, thus 
maintaining an (approximately) uniform mesh in the BZ. 
The ground state of bulk GaAs, which is also considered 
for reference, is computed using a 12x12x12 sampling. 

We compare the density of states at the KS level^, 
for the four systems considered, in figure [2j All the 
Mn-containing systems have a metallic character, as ex¬ 
pected. The case of undoped bulk GaAs is also shown 
for comparison. Doped systems clearly display a nonzero 
DOS at the Fermi level, with a characteristic Mn-related 
structure lying in the top region of the bulk bandgap. 
The 50% doped structure (panel (a)) clearly deviates 
from the behaviors shown by lower-doping systems, both 
at the level of the DOS and in term of Fermi surface 
topology. This difference can be clearly seen in the band 
structure, where the number of bands appearing to cross 
the Fermi level increases from three (see e.g. Fig. ED, 
for the lower doping cases, to four (see Fig. S]) in the 
50% doped system. At the same time, the DOS of the 
50% doped case displays a clear reduction of the spectral 
weight around 2.5 eV above E^?, typical of bulk GaAs, 
in favor of lower energy structures. 

Moreover, the total magnetic moment per unit cell is 
an integer multiple of jub (namely, 4//^) in all the struc¬ 
tures except the 50% doped one, where it drops slightly 
below the value of 4 A non-integer value of fi / /a b 
implies that the system cannot be an half metal, as also 
appearing from spin projected DOS (Eig. [2]), where both 
the spin-up and spin-down component of the correspond¬ 


ing DOS are non-negligible at the Eermi energy. In the 
inset of each panel we display the spin-down component 
of the DOS, zoomed in the vertical scale. 

The 12%, 17%, and 25% DEHs present a gap in the 
spin-down projected DOS, which is zero at the Eermi 
energy. At variance, the spin-up DOS of all these DEHs is 
clearly nonzero at the Eermi energy. Hence, the electron 
population at the Eermi energy in the 12%, 17%, and 
25% DEH is fully spin-up polarized, and these systems 
are half-metals, similarly to the case of uniformly Mn 
doped GaAs compounds^. 

In view of the qualitative difference of the 50% doped 
system, in the following discussion and in optical proper¬ 
ties calculations only materials with integer ji / jiB will 
be considered, and compared with the uniformly doped 
systems. 

A more detailed analysis of the band character can be 
performed by decomposing the DOS via a projection on 
the different atomic species and angular momenta. The 
major contribution to the metallic bands comes from Mn 
d orbitals and As and Ga p orbitals. Eig. [5] shows the As 
p (top panel), Mn d (middle panel), and Ga p (bottom 
panel) contribution for the three systems considered (i.e. 
three, five and seven layers super-cell). The As and Ga p 
contribution to the undoped GaAs DOS is also shown for 
comparison (black line). The p As and Ga contributions 
are also decomposed according to the distance of the con¬ 
sidered atom from the Mn layer, with the thick contin¬ 
uous line corresponding to the As (Ga) nearest to Mn. 
As shown in figure [5l the presence of Mn mainly affects 
the first two closest layers of As (labeled with 1 and 2 in 
Eig. [TJ thick continuous and dashed lines in Eig. [5]) and 
the first closest layer of Ga (same conventions). They 
all display a peak above the Eermi energy at the same 
energy of the d orbital of Mn. These peaks indicate an 
hybridization of d orbitals of Mn with p orbitals of the 
As and Ga atoms, or better with the sp^ backbone of 
the GaAs lattice. Indeed the Mn atoms, replacing Ga 
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Figure 5. (color online) Projected density of states of the 
studied MnGaAs heterostructures. The projections on As p 
orbitals (panels al-a3), Mn d orbitals (panels bl-b3) and Ga 
p orbitals (panels cl-c3) are shown for the structures with 3 
(blue), 5 (green) and 7 (red) GaAs layers between the Mn 
layers, respectively. The projection on As and Ga species is 
atom resolved. Atoms are selected according to their distance 
from the Mn layer, as 1st nearest neighbor (thick continuous 
line), 2nd (thick dashed line), 3rd (dotted line) and 4th (thin 
dashed line). The case of bulk GaAs is also shown for com¬ 
parison (thin continuous black line). 


atoms, are “forced” to create an sp^ hybridization. How¬ 
ever, while Ga has three electrons in the s plus p orbitals, 
Mn only has two. There is then a competition between 
filling the hole in the sp^ backbone and keeping five va¬ 
lence Mn electrons in its d orbitals. This competition is 
represented as two resonant configurations in Fig.[6l The 
dashed red arrow indicates the electron shared between 
sp^ and d orbitals. In both configuration the magnetic 
moment is 

Remaining As and Ga atoms show a PD OS contribu¬ 
tion very similar to the one in bulk GaAs^ confirming 
that they are not substantially influenced by the pres¬ 
ence of the Mn dopant. This result confirms that the 
ferro-magnetic interaction between the Mn layers be¬ 
comes much weaker if they are separated by more than 
3 GaAs layers. However also at greater distance a small 
projection of the state immediately above the Fermi level 
remains on the As and Ga orbitals, suggesting that a 
smaller ferro-magnetic interaction is still present also in 
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Figure 6. Schematic electronic diagram for the As-Mn bind¬ 
ing represented as two resonant configurations. Each Mn 
orbital has four first nearest neighbor As atoms. In first ap¬ 
proximation, an hole is shared between the four Mn-As sp^ 
bonds orbitals and the Mn d orbitals. In practice the hole 
delocalize also on the other sp^ orbitals of the doped GaAs 
and is responsible of the ferro-magnetic interaction between 
different Mn layers. 


the 12% doped system, where the Mn layers are sepa¬ 
rated by 7 GaAs layers. 


III. OPTICAL PROPERTIES OF Mn DOPED 
GaAs 


Optical absorption is proportional to the imaginary 
part of the macroscopic dielectric function, cm (<^)- The 
latter can be obtained from microscopic quantities such 
as the response function XG,G'iQ^ Response func¬ 

tions can in turn be computed on the basis of transition 
matrix elements between electronic states (both inter¬ 
band and intraband in the case of metals). Several levels 
of theory can be adopted, ranging from the simplest ap¬ 
proach, i.e. the bare independent particle-random phase 
approximation (IP-RPA), to more complete approaches 
where, e.g., local fields and/or excitonic effects are taken 
into account. In the long wavelength limit, which is ap¬ 
propriate for UV-VIS spectra, cm can be expressed as: 


cmH = 1 - lim(1) 

Involving only the G = 0, G' = 0 “head” of the appro¬ 
priate response function XG,G'iQ^^) written in reciprocal 
spacei^. 

In this work, we use XcG^i^^ response func¬ 

tion constructed including the local-field (LF) effects: 










































































































5 


XGG^(95w) —Xgg>(9,<^) + 

+ XGGl(9,‘^)fG^G2(9,‘^)XG2G'i9^‘^) 

G1G2 

( 2 ) 

Here /q^q^ (quj) = 6vh[p]/Sp, the functional derivative 
of the Hartree potential, is the bare Coulomb interaction 
and corresponds to an electron-hole exchange terrr^. 

Neglecting /" one has xgg'(<?,w) = [x^^]gg'( 9 , w), 
i.e., the response function in the IP-RPA approximation. 

The quasi-particle (QP) response function an be ex¬ 
pressed in terms of single-particle electronic eigenstates 

eigen-energies as: 

XgG'(9,^) = -^ (/nfe - fmk') 

nk mk' 





^QP obtained from KS eigenvalues applying a gap¬ 
opening/ band-stretching correction which simulates the 

self-energy effects in semiconductors/metals, and 

are taken as unperturbed KS eigenstates as 

plemented in the YAMB02^ code. 

Here we report the optical properties of our four, six 
and eight layers SC of Mn doped GaAs. The absorption 
spectra of the undoped system is also shown for reference 
as for the case of the DOS. To obtain converged spectra a 
k-points sampling of the BZ with 24 x 24 x n mesh is used, 
with n = 8,6,4 for the four, the six layers and the eight 
layers respectively, considering nv valence bands and nc 
conduction bands in total. The Drude term is modeled 
with a plasma frequency ujp = 2.75eK and relaxation fre¬ 
quency of 0.7eK and the QP corrections are included as 
a 10% stretching of the band structure, as in Ref. [I 3 . 
LF effects are found to converge including G vectors in 
Eq.[2]up to a kinetic energy a cutoff of ARy. The absorp¬ 
tion of bulk GaAs is instead computed with a double-grid 
sampling of the BZ2^ with a regular 12x12x12 mesh 
shifted from F plus a random mesh of 2000 k-points. In 
this case the well established valued! of 0.8 eV opening 
of the band gap has been used for the QP correction. 

In Figs. [TllH] the absorption spectra, for the 

12, 17 and 25% doped materials is shown. The case of 
bulk (undoped) GaAs is in panel (a) of both figures. 
Dashed lines show the IP-RPA spectra for comparison. 
The RPA spectra including LF effects differ from IP-RPA 
ones because the intensity of the absorption is reduced. 
The intensity reduction is peak-specific, and is due to the 
well known depolarization effect, which is usually more 
important in non uniform systems. However, LF do not 



Figure 7. (color online) Computed UV-VIS absorption spec¬ 
tra for light polarization parallel to the Mn layers. Panel (a) 
shows the bulk GaAs case as reference, and Mn (5-doped het¬ 
erostructures with 12.5%, 16.7%, and 25% Mn concentration 
are shown in panels (6), (c), and (d), respectively. The dashed 
lines and continuous lines represent the IP-RPA and RPA re¬ 
sults, respectively. Theoretical spectra are compared with 
experimental data for bulk GaAs (circles in panel (a), from 
Ref. [23]) and for a uniformly Mn-doped GaAs sample at low 
doping concentrations (dots in panel (6), from Ref. [23]). We 
also compare with the FLAPW ab-initio results of Ref. 0, 
computed for a uniformly doped system with a concentration 
of 6.25% (dot-dashed line). 


change substantially the position of the peaks. 

We first consider incident light with polarization on 
the layers plane (Fig. [7j). The intensity of the main E2 
peak around 5 eV is lower in DEHs than in bulk GaAs, 
and it decreases by increasing the Mn concentration. In 
fact this peak becomes a double peak by effect of the Mn 
d-doping (panel (5)), and more structures appear by de¬ 
creasing the distance between Mn layers (panels (6), (c) 
and (d)). We compare in panel (5) our IP-RPA results for 
the 12% delta-doped systems (seven layers of Ga between 
the Mn layers) with experimental data (circles in panel 
(6)) for an uniformly doped Gai-xMuxAs sample with 
a; = 0.066 (Ref. [1^). The theoretical underestimation of 
the intensity of the 5 eV (E2) peak is likely due to the 
different doping concentration between theory and ex¬ 
periment. Indeed the intensity of the E2 peak decreases 
increasing the dopant concentration, both experimentally 
fRefs . [^ andl^f and theoretically (present work, and 
Ref. [lO])- We notice that the double peak structure at 
^ 5 eV of Eig.[3is peculiar of the d-doped GaAs absorp¬ 
tion spectra in polar geometry, since it is absent in the 
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Figure 8. (color online) Computed UV-VIS absorption spec¬ 
tra of bulk GaAs (panel (a)) and Mn (5-doped heterostruc¬ 
tures with 12.5% (panel (5)), 16.7% (panel (c)), and 25% 
(panel (d)) Mn doping, for light polarization perpendicular 
to the Mn layers (grazing incidence). As in figure [71 the 
dashed lines and continuous lines represent IP-RPA and RPA 
results, respectively. 


computed spectra with light-polarization perpendicular 
to Mn layer (Fig. [8]), as well as in the experimental data 
for randomly distributed Mn doping. The double peak 
also appears when the absorption spectra is computed 
without SO as shown in Fig He. 2. 

The influence of local field effects in the case of grazing 
incidence (Fig. |8j) is more pronounced for energies above 
5 eR, where LF appear to suppress the peak intensity 
more than at lower energies. 

On the other hand, in both Figs. 0 and [8] the inten¬ 
sity of the theoretical ^ 3 eV is underestimated, due to 
our neglect of excitonic effect The latter, however, 
are expected to be strongly suppressed in systems with 
efficient free-carrier screening, such as in metals or even 
half-metals, with respect to the case of a pure semicon¬ 
ductor. In fact the peak at 3 eV in the experimental data 
for uniform Mn doping (dots in panel (6) of Fig. [7] ) is 
clearly less prominent than in bulk GaAs. 

The overestimation of cm at very low energy in our 
simulation can be attributed, at least in part, to an 
overestimation of intraband transitions, described by the 
Drude term. The discrepancies between experimental 
data and our results can also be partially accounted for 
by the different type of doping (uniform vs ^-doping). 

Moreover experiments are performed at diluted Mn 
concentration (few percents or less). When the concen¬ 
tration of Mn in uniformly doped samples is sufficiently 


low, the overlap of impurity wave-functions is negligible 
and the compound is semiconducting. At variance, in 
^-doped samples, the conductivity in the Mn layers re¬ 
mains metallic also at overall low doping concentrations. 


A. Kerr effect and reflection magneto—optical 
dichroism 


MOKE experiments evaluate the difference in optical 
reflectivity for left and right circularly polarized light. In 
practice, since a linearly polarized light can be thought as 
the result of the superposition of left and right circularly 
polarized waves having the same intensity and phase, a 
nonzero difference in the dielectric response of the ma¬ 
terial can be detected as a change in the polarization 
properties of the reflected light. 

In this way, Kerr parameters are defined in terms of 
the Kerr rotation Ok and the Kerr ellipticity rjK , the first 
describing the angle between the major axis of the ellipse 
with respect to the original linear polarization direction, 
and the second one being related to the ratio between 
the two ellipse axis. 

In the present manuscript we consider the MOKE in 
the so called polar geometry, where the light propagates 
along the direction of the magnetization of the sample, 
which we identify here as z, perpendicular to the sample 
surface. The polarization hence lays in the xy plane. We 
now introduce the complex refractive index n± = yTT, 
where (e_) is the dielectric function for the right (left) 
circularly polarized light. Then following Ref. [gI the Kerr 
parameters are defined as: 


Ok {^) + 'iVK (<^) 


. nj^iuj) — n-iuj) 

% - . 

nj^{uj)n-{uj) — 1 


( 4 ) 


e± can be derived, for a cubic system, from the xx and 
xy elements of the dielectric tensor, according to the ex¬ 
pression: e± = exx The magneto-optical response 

can hence be computed starting from the diagonal and 
off-diagonal components of the dielectric tensor. Indeed, 
if exy <C €xx^ which is the common case, EqH] is well 
approximated by: 


Ok (^) + ir]K (^) 


_ ^xy (^)_ 

[^XX (^) 1 ] \/ ^xx (^) 


( 5 ) 


which holds for small Kerr angles. 

Remarkably, in such differential measurements, the 
contribution stemming from the off-diagonal elements of 
the dielectric tensor, which is usually negligible in bare 
absorption or reflectivity experiments, becomes impor¬ 
tant. Cxy can be obtained from the response function 
by generalizing Eq. [H In practice, we evalu¬ 
ate the dielectric tensor by means of the YAMBO^ code, 
using the same approach described in Ref. [TsI, where 
€a/ 3 , is computed, at the RPA-IP level, as 











Kerr angle 


Ellipticity 
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47re^ 

eMa/3 = 1 - lim — (QaQisi^) ■ (6) 


exy must be converged against the same parameters 
needed to converge the absorption. Convergence in the 
sampling of the BZ requires a 24 x 24 x n grid where the 
number of k-points along the 2; direction are n = 12, 8, 6 
for the 4, 6 and 8 layers super-cell respectively i.e. a 
slightly finer sampling along the 2; direction, if compared 
to the one used for absorption. The reason stems from 
the fact that exy is two orders of magnitude smaller than 
exx> The knowledge of also allows one to compute 
the reflectance R±{uj) at normal incidence, defined as 
the square modulus of the complex reflectivity, both for 
right (+) and left (—) circularly polarized light^: 


n±{uj) — 1 
n±{uj) + 1 


(7) 


Following Ref. [H, the latter quantities determine the 
reflectance magnetic circular dichroism (R-MCD) spec¬ 
trum as: 


R 


MCD 


M = 


I R+-R- 

2 R-\- T R— 


(8) 


which is more easily measured experimentally than the 
MOKE parameters. In the limit of small Kerr angles one 
can prove, after some tedious but straight toward alge¬ 
bra, that Rmcd{^) and Kerr ellipticity are brought to 
coincide: 


Rmcd{^) ~ (9) 

In Fig.[9]we plot the computed Kerr parameters Ok and 
vjK for the considered DFH. We also show, by compar¬ 
ison, some theoretical and experimental literature data 
for uniformely doped MnGaAs. 

We first observe that there is no clear dependence of 
the computed Kerr parameters on the Mn doping con¬ 
centration. This finding is in agreement with previous 
calculations for the uniform doping case^i^. Our calcu¬ 
lations shows that, as it happens for the absorption, the 
computed MOKE is similar to the one obtained for the 
uniformly doped case. In Fig. [9] (bottom panels) the Kerr 
rotation is compared with all-electrons calculations for 
the uniform system at the LDA leveUS. Despite this dif¬ 
ference in the calculations, different structures (uniform 
doping and DFH) and different doping concentrations all 
calculations and experimental data shows at least one 
common dominant feature, that is a negative Kerr ro¬ 
tation at low energy. In the experimental data for uni¬ 
form dopin g^^i^^ this negative peak in the Kerr rotation 
is blue-shifted by 0.6 eV with respect to our results 
(see Fig.O top left panel). This shift also appears in the 
computed Kerr ellipticity/R-MCD data (top right panel) 
and in DFT calculations for the uniform doping case. It 



Figure 9. Kerr parameters computed at the IP-RPA level 
for three Mn (5-doped GaAs heterostructures. Central pan¬ 
els show the spectrum of the 12.5% Mn DFH (red line); the 
results for the 16.7% and 25% Mn concentrations are shown 
in bottom panels (green and dashed blue lines, respectively). 
In the top panel, we show, for comparison, several experi¬ 
mental data for uniformely doped samples: black dots are 
for a 6% Mn sample considered in Ref. [s^, the red dots 
are for a 3% Mn case (Ref. 0) , and violet dots refer to a 
2% Mn sample (Ref. [2^) Finally blue empty dots in right 
top panel are experimental data for the Reflection MCD in a 
uniformely doped sample with 3% Mn from Ref. [ 2 ^ . In bot¬ 
tom panels, we report as dashed lines the results of FLAPW 
(LDA-hstretched bands) calculations for an uniformly doped 
systems, from Ref. 0 


maybe due to inadequacy of the 10% band stretching 
approach to mimic QP corrections, in particular for the 
Mn-d empty states. The latter, being only 1 eV above 
the Fermi level, are essentially unaffected by the stretch¬ 
ing approach. 

We can thus conclude that at least the first peak in the 
Kerr rotation remains substantially unchanged for dop¬ 
ing concentrations varying from 12.5% to 25%, and it is 
barely affected by the details of the geometrical distribu¬ 
tion of the Mn atoms. 
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Figure 10. Real and imaginary part of the diagonal and 
off-diagonal components of e for the Mn (5-doped GaAs het¬ 
erostructure corresponding to a 12.5% doping (Ref. HI) 


As it appears from Eq. [5l a peak in the Kerr rotation 
can either have an “optical” origin, i.e. it stems from 
a vanishing value of the denominator, or a “magneto¬ 
optical” one, i.e. it arises from maximum/minimum of 
exy{oo) in the numerator. For the uniformly doped system 
the minimum close to 1 eV has been reported to be of 
“optical” origir^ii^. 

To analyze if this is the case also for the DFH, we plot 
in Fig. [To] the diagonal and the off-diagonal dielectric 
function. We clearly see that Im[exy] has a maximum at 
1 eV and also that both Re[exx] and Im[exx] are small 
in the same region. So the corresponding peaks of the 
Kerr rotation, in the present case, have both optical and 
magneto-optical origin. 

Comparing the computed Kerr parameters with the 
computed DOS and its projections (see Figs. 11151 ) , we see 
that a key role is likely played by the empty state induced 
by the Mn doping, as discussed in sec. [11] (see also Fig. [6}. 
Indeed we note that this state is located about 1 eV above 
the Fermi level. Thus the maximum (in absolute value) of 
the rotation is likely due to transitions from the valence 
band to this state, i.e. it is due to a transition involving 
both d and sp^ states. 

We thus suggest that in DFH the MOKE could be 
substantially affected by electron/hole doping. Indeed, 
varying the holes concentration, the semi-metallic char¬ 
acter of the system can be substantially affected, leading 
to a very different absorption spectra, in particular in the 
region close or below 1 eV as already suggested for the 
uniformly doped systen>^. This conclusion is also sup¬ 
ported by previous all-electrons results for the uniform 
doping case, where the Kerr spectrum is shown to change 
significantly when As anti-site defects are considered to¬ 
gether with Mn doping^. 


Faraday angle 


Ellipticity 



[eV] 


Figure 11. Faraday parameters / transmission MOD com¬ 
puted at the IP-RPA level for three Mn (5-doped GaAs het¬ 
erostructures, corresponding to 12.5%, 16.7%, and 25% Mn 
concentrations. Black empty dots are experimental data for 
uniformly doped samples with 1% Mn, from Ref. [s^. 


B. Faraday effect and transmission 
magneto—optical dichroism 


While the Kerr effect and Reflection MCD imply mea¬ 
suring the light reflected by a magnetic material, the cor¬ 
responding effect for light transmitted through a mag¬ 
netic medium is known as Faraday effect. In polar ge¬ 
ometry (the one considered in the present study), when 
a linear polarized light propagates parallel to the magne¬ 
tization of a ferromagnet, a Faraday rotation is observed 
as a rotation of the polarization plane of the transmit¬ 
ted light {Of)- In term of circularly polarized reflective 
indexes, the Faraday rotation Op and the Faraday ellip¬ 
ticity r]F read^i 


Of {uj) -h ir]F {uj) = (n+(cj) - n_(cj)) L (10) 

where L is the sample thickness and c the speed of light. 
We evaluated the Faraday parameters per unit length for 
the Mn (5-doped GaAs in a similar way as done for Kerr 
parameters; our results are displayed in Fig. [TTl As the 
Kerr parameters are related to the R-MCD, the Faraday 
parameters are related to the transmission MCI3^ (T- 
MCD): 


Tmcd{^) = F^FpF^ ( 11 ) 
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which, at normal incidence, is proportional to A(a(cc;), 
i.e. the difference of the absorption coefficient 
= ‘2ujlm[n±{uj)]/c between the right and the left 
circularly polarized light: 

A(a(co’) = (a+(co’) — a-{uj) = 

1j 

Thus we can compare the results of our simulations with 
T-MCD experimental data for diluted concentration^ in 
randomly Mn-doped GaAs. 

First thing that we observe is that, despite the physics 
of the Faraday effect / T-MCD is very similar to the on 
of the Kerr effect / R-MCD, in this case the results are 
slightly more sensitive to the doping concentration. At 
least the region above 2 eV is affected going from 12.5% 
to 25% doping. However, as in the case of the Kerr, there 
is a structure slightly above 1 eK, which is substantially 
unchanged varying the doping concentration. 

To investigate possible applications of DFH as Faraday 
rotators, in the following we compare some of our com¬ 
puted optical and magneto-optical data with the corre¬ 
sponding quantities of bulk transition metals, which are 
typical textbook examples for the Faraday effect 

We consider the 4-layers DFH and the peak at ^ 1 
eV: according to our simulation, the Faraday rotation is 
, Op ^ significantly larger than the typical val¬ 
ues for Fe, Co, and Ni, (35, 36, 10 respectively)^ 

(taking into account the 0.6 eK red shift of our computed 
spectra with respect to the experimental ones, these val¬ 
ues are taken at A = 820 nm, i.e. at about 1.5 eV). 
On the other hand, DFH have a lower optical absorption 
coefficient, which is desirable for efficient transmission. 
The optical penetration depth, defined as the length over 
which an electromagnetic wave of given frequency en¬ 
counters a 1 /e reduction of its intensity -proportional to 
the inverse of the absorption coefficient- is about 36 nm 
for our 4 layers DFH. For comparison, the optical pene¬ 
tration depth of transition metals are of the order of 15 
nm. 


IV. CONCLUSIONS 

In conclusion, our results for Mn (5-doped GaAs show 
that such digital ferromagnetic heterostructures have op¬ 
tical absorption and magneto-optical spectra displaying 
both similarities and differences with respect to previ¬ 


ously studied uniformly doped MnGaAs systems. The 
50% Mn doping case (one layer of GaAs between MnAs 
layers) corresponds to a system with qualitatively differ¬ 
ent properties, having lost the half-metallicity. All other 
structures, corresponding to 12-25% Mn doping, have op¬ 
tical absorption spectra which change gradually with the 
Mn concentration. The resulting spectra compare well 
with the available experimental data for the low-doping 
(< 12%) system and theoretical data for the higher dop¬ 
ing (25%), for the uniformly doped case. With increasing 
Mn concentration the typical bulk GaAs spectral shape 
becomes less and less structured. In agreement with the 
experimental findings the peak at 3 eV is strongly sup¬ 
pressed, likely because of the hole doping induced by the 
Mn impurities, which screens the electron-hole interac¬ 
tion responsible for the creation of the excitonic peak. 
The peak at 5 eV is instead less affected. 

Although the absorption spectra of the (5-doped mate¬ 
rials share the essential features of those of uniformly Mn 
doped GaAs at similar concentrations, the effect of local 
fields is enhanced, due to the stronger non-uniformity of 
the ^-doped system. This effect is particularly evident 
when the light polarization has a large component along 
the growth direction, while it is less important when the 
polarization lays in the plane of the Mn layers. MOKE 
spectra do not show a clear trend with the doping con¬ 
centration: essentially, the maximum value of Kerr rota¬ 
tion is constant, and quite similar to that of the uniform 
doped material^i^. Although the Faraday ellipticity is 
more sensitive to the doping concentration, also in this 
case the low energy spectrum remains substantially un¬ 
changed moving from the high doping case (25%) to lower 
doping. On the other hand, the important role of the hole 
shared between the Mn d and GaAs sp^ states suggests 
that the studied DFH could be very sensitive to a further 
doping with electrons. - 
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